A New Giant Stellar Structure In the Outer Halo of M31 
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^ , The Sloan Digital Sky Survey has revealed an overdensity of luminous red 

giant stars ~ 3° (40 projected kpc) to the northeast of M31, which we have called 
\ Andromeda NE. The line-of-sight distance to Andromeda NE is within ~ 50 kpc 

of M31; Andromeda NE is not a physically unrelated projection. Andromeda 
Q\ ', NE has a g-band absolute magnitude of ~ —11.6 and central surface brightness 

of ~ 29 magarcsec -2 , making it nearly two orders of magnitude more diffuse 
than any known Local Group dwarf galaxy at that luminosity. Based on its 
q ! distance and morphology, Andromeda NE is likely undergoing tidal disruption. 

\ Andromeda NE's red giant branch color is unlike that of M31's present-day outer 

P-T disk or the stellar stream reported by Ibata et al. (2001), arguing against a direct 

link between Andromeda NE and these structures. However, Andromeda NE has 
a red giant branch color similar to that of the Gl clump; it is possible that these 
structures are both material torn off of M31's disk in the distant past, or that 
. ' these are both part of one ancient stellar stream. 

Subject headings: galaxies: general — galaxies: evolution — galaxies: stellar 
content — galaxies: dwarf — galaxies: individual (M31) — Local Group 



1 Max-Planck-Institut fur Astronomie, Konigstuhl 17, D-69117 Heidelberg, Germany; zucker@mpia.de 

2 Astronomisches Institut, Universitat Basel, Venusstrasse 7, CH-4102 Binningen, Switzerland 

3 Astronomy Department, University of Washington, Box 351580, Seattle WA 98195-1580 

department of Astronomy, New Mexico State University, 1320 Frenger Mall, Las Cruces NM 88003-8001 

5 Princeton University Observatory Princeton, NJ 08544 

6 Apache Point Observatory, P.O. Box 59, Sunspot, NM 88349 

7 Astronomy and Astrophysics Department, University of Chicago, 5640 S. Ellis Ave, Chicago, IL 60637 



- 2 - 



1. Introduction 

The study of dwarf galaxies offers special insight into the physics of galaxy formation 
and evolution. There are many fewer dwarf galaxies than predicted by a direct scaling of 
the dark matter halo mass function (e.g., White & Frenk 1991). A number of theoretical 
explanations for this have been explored: feedback from supernovae leading to disruption of 
dwarf galaxies (e.g., Dekel & Silk 1986) or strong reduction of their luminous/dark matter 
ratio (e.g., Hayashi et al. 2003), reduced substructure in warm dark matter haloes (e.g., 
Moore et al. 2000), enhanced satellite disruption and reionization (e.g., Bullock, Kravtsov, 
& Weinberg 2001), or suppression of dwarf galaxy formation through photoionization (e.g., 
Somerville 2002). In this context, the dwarf galaxy content of the Local Group plays a 
special role. The luminosity function can be probed down to very faint limits, where the 
discrepancies with the models are largest (Moore et al. 1999). Furthermore, the detailed 
star formation histories of the dwarfs can be estimated through analysis of the colors and 
magnitudes of their resolved stellar populations. Coupled with the exquisitely deep surface 
brightness limits achievable through analysis of resolved stellar populations, one can explore 
the physical processes which shape dwarf galaxy evolution directly, such as ram pressure 
stripping (e.g., Grebel, Gallagher, & Harbeck 2003) or tidal disruption by the giant galaxies 
in the Local Group (e.g., Ibata, Gilmore & Irwin 1994). 

Nevertheless, constructing a complete census of faint stellar structures (dwarf galaxies, 
disrupted dwarfs, stellar streams) has proved challenging, owing to the twin constraints 
of imaging depth and large area coverage (see, e.g., Armandroff, Jacoby, & Davies 1999; 
Ferguson et al. 2002). In this letter, we report the discovery of a faint, extended stellar 
feature in the outer halo of M31, using a special scan from the Sloan Digital Sky Survey 
which reaches depths sufficient to detect stars at the tip of the red giant branch (TRGB). 
We use these data to estimate the feature's physical properties and to try to determine its 
true nature. We assume a distance to M31 of 760 kpc (van den Bergh 1999). 

2. Observations and Data Analysis 

The Sloan Digital Sky Survey (SDSS) (York et al. 2000) is an imaging and spectroscopic 
survey that will eventually cover about 1/4 of the sky. Drift-scan imaging from the SDSS 
2.5m telescope in the five SDSS bandpasses (u,g,r,i,z) (Fukugita et al. 1996; Gunn et al. 
1998; Hogg et al. 2001) is processed through data reduction pipelines to measure photometric 
and astrometric properties (Lupton et al. 2002; Stoughton et al. 2002; Smith et al. 2002; Pier 
et al. 2003) and to identify targets for spectroscopy. M31 was observed on 5 October 2002, 
during time when primary survey areas were inaccessible. The data were obtained in two 
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overlapping strips oriented along the major axis of M31, covering an area of ~ 18° x 2?5. 
Conditions were photometric, but the seeing varied between 1" — 2" FWHM in r-band. The 
raw data were processed and catalogs constructed using the same pipeline as Data Release 
1 (DR1) (Abazajian et al. 2003). We selected only objects within the coordinate range 
assigned to each strip, thus eliminating duplicate object detections. In the following, all 
references to dereddening and conversion from SDSS magnitudes to V, I magnitudes (for 
literature comparison) utilize Schlegel, Finkbeiner, & Davis (1998) and Smith et al. (2002), 
respectively. 

The spatial density of SDSS-detected stars with 21 < i < 22 is shown in Fig. 1. The 
central regions of M31, its disk, and its bright satellites M32 and NGC205 are too crowded 
for the standard SDSS photometric pipeline. The slight horizontal "striping" is the result 
of seeing variations; under better seeing conditions, more stars were detected. This effect is 
most pronounced at faint magnitudes. Complex stellar substructure is evident in the outer 
disk and halo of M31 (Fig. la). Essentially all of the morphological features visible in 
Figures 2 and 3 of Ferguson et al. (2002) can also be seen in the SDSS data, although the 
giant stellar stream to the southeast (Ibata et al. 2001; Ferguson et al. 2002; McConnachie 
et al. 2003) is truncated by the limited width (2?5) of the SDSS scan. 

Adjusting the grayscale in Fig. la to enhance faint details, we obtain Fig. lb, in which 
an overdensity around (X,Y) ~ (0.5,3.1), or 00 h 52™0 + 44°06' (J2000), becomes apparent. 
The seeing, rather fortuitously, was almost constant in the vicinity of this overdensity, min- 
imizing the striping in this area. The stellar overdensity has its highest concentration at a 
angular separation of ~ 3?3 from the center of M31, and although it is adjacent to the area 
mapped by Ferguson et al. (2002), it is not covered by the data shown in that work. 1 

In order to place constraints on the nature and distance of this new stellar overdensity, 
we dereddened the data, transformed g, r, i photometry to V, I magnitudes, and constructed 
color-magnitude diagrams (CMDs) of a region in the outskirts of M31, of the new feature, and 
of a nearby control field (Fig. lc). The left side of Fig. 2 shows the Hess diagrams of the M31 
periphery (between the two ellipses in Fig. 1), the control field, and the foreground-corrected 
M31 Hess diagram. The right side of Fig. 2 shows the same for the new, diffuse stellar feature. 
The new stellar feature has a Hess diagram which is qualitatively similar to that of M31, 
placing it roughly at the distance of M31, and firmly establishing its extragalactic nature. 



1 Subsequent to initial submission of this work, Lewis et al. (2004) presented a map which clearly shows 
this feature, although the authors did not comment on its significance. 
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3. A New Stellar Structure: Andromeda NE 

The stellar overdensity revealed in Fig. lb, which we will call Andromeda NE (And 
NE) in what follows, covers more than half a square degree. To minimize foreground star 
contamination in such an extended, faint stellar feature, we selected stars in those Hess 
diagram bins from the bottom left panel of Fig. 2 with S/N > 3.0 (i.e., those likely to 
be at the distance of M31) for further analysis (cf. Odenkirchen et al. 2001). The SDSS 
data only include stars near the TRGB (2.6 arcmin -2 at the densest part of the feature), 
and so to estimate And NE's peak surface brightness and total luminosity we carefully 
accounted for this effect. Following the method of Kniazev et al. (2004), we measured the 
total luminosity of all selected stars in And NE in circular apertures, minus an average 
stellar background level determined from nearby external fields. Analyzing SDSS data for 
the Pegasus dwarf irregular (DD0216), a galaxy at a distance comparable to that of M31, we 
fitted surface brightness profiles to both stars and to the total integrated light in its outer 
regions, determining offsets of ~ 2.8, 2.5, and 2.0 mag in g,r, and i between the 'resolved- 
source' and total surface brightnesses (an i-band correction of 1.8 mag was derived from the 
i-band stellar luminosity function of the Draco dwarf spheroidal given in Odenkirchen et al. 
2001, confirming that these estimates are reasonable for a range of likely stellar populations). 
We then applied these corrections to the And NE measurements, obtaining central surface 
brightnesses and total magnitudes of 29 ± 0.4, 28.3 ± 0.3, and 27.2 ± 0.3 mag arcsec -2 , and 
12™8, 11™4, and 11™0 in g, r, and i-band, respectively. If we assume that And NE is at 
approximately the same distance as M31 (justified in the next paragraph), and correct for a 
median foreground extinction of A g ~ 0.32 (A v ~ 0.28), it has M g ~ -11.9 (M v ~ -12.7); 
Local Group dwarfs of comparable luminosity have V-band surface brightnesses nearly two 
orders of magnitude higher than And NE (see, e.g., Grebel, Gallagher, & Harbeck 2003). It 
is clear that And NE, if it is a galaxy, is a highly unusual one. 

We estimated the distance to And NE from the position of the TRGB (cf. Lee, Freedman, 
& Madore 1993) in And NE. We generated cumulative dereddened /-band luminosity func- 
tions for stars within the two circular regions shown in Fig. lc, and subtracted the control 
field luminosity function from that of And NE. The field-star-subtracted luminosity function 
rises noticeably at / ~ 20.4; for RGB metallicities in the range —2.2 < [Fe/H] < —0.7, 
Mj^trgb ~ —4.0 ±0.1 (Da Costa & Armandroff 1990; Salaris & Cassisi 1998), yielding a 
distance modulus of (m — M) ~ 24.4 ± 0.15 (760^50 kpc) for And NE. Thus we conclude 
that And NE is at approximately the same distance as M31. However, the precision of this 
result is limited by the small number of stars at the TRGB in And NE and the possibility 
that And NE is extended along the line of sight; significant improvement will likely only be 
possible using, e.g., more numerous red clump or horizontal branch stars. 
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What is And NE? To give more insight into its nature, we compared the Hess diagrams 
for the red giant branch (RGB) in And NE and in three other prominent stellar features: 
the giant stream, the northern spur, and the Gl clump. For each field, we dereddened the 
SDSS data, transformed g,r, % stellar magnitudes to V and J, subtracted an appropriately 
scaled control field Hess diagram, and finally divided by the square root of the sum of the 
target and control field Hess diagrams to estimate S/N in each bin. Fig. 3 shows the 
results, with overplotted RGB fiducials for Galactic globular clusters spanning a wide range 
of metallicities. And NE has a "blue" RGB (following the nomenclature of Ferguson et al. 

2002) , unlike the "red" RGBs shown by the northern spur and the giant stream, from which 
we infer that And NE has a lower mean metallicity. For this reason, we believe it unlikely 
that And NE is an extension of either the nearby northern spur or the giant stream. 

Interestingly, And NE and the Gl clump have very similar RGB morphologies, despite 
the large differences in their stellar densities and projected positions (70 kpc on opposite 
sides of M31's disk). Moreover, And NE appears to have a faint extension towards the disk 
of M31. Taken together, these data suggest two possibilities for And NE. The first is that the 
Gl clump and And NE are both material torn from M31's disk by an ancient interaction; this 
possibility is supported by the disk-like kinematics of RGB stars in the Gl clump (Reitzel et 
al. 2004), but is argued against by the strong difference in RGB morphology between M31 
disk stars and And NE (M31's disk appears more metal-enriched than And NE or the Gl 
clump; see Figs. 2 and 3). Secondly, it is possible that And NE and the Gl clump are both 
projections along the orbit of an ancient, diffuse stellar stream roughly in the plane of M31's 
disk. The morphology of both features and the lack of similarity between disk and And 
NE stellar populations support this possibility. The disk-like kinematics of stars in the Gl 
clump would appear to argue against this scenario, but such a coincidence is not impossible; 
the Monoceros stream around the disk of the Milky Way (Newberg et al. 2002; Yanny et al. 

2003) exhibits prograde, disk-like kinematics (e.g., Penarrubia et al., in preparation). 

Having discussed possible stellar stream and stellar debris scenarios for And NE, we turn 
to the issue of whether And NE is a bound, self-gravitating entity (i.e., a new satellite). We 
can calculate what mass-to-light (M/L) ratio would be required for it to be gravitationally 
bound by using an approximation for tidal disruption for a cluster of particles in a circular 
orbit: 3M M3l /Dl l31 > Mamne/ -RindNE; where M M3 i and M AndNE are the total masses of 
M31 and And NE, -Dm3i is the distance of And NE from M31, and i?AndNE is the radius 
of And NE. The masses are unknown; however, we can calculate the total M/L of And NE 
that would be required to bind And NE, in units of M31's M/L at radii > 40 kpc. The 
integrated magnitude of And NE is Vamne ~ 12.0, and using Vm3i = 3.44 from the RC3 
catalog, and assuming that M31 and And NE have the same reddening and lie at the same 
distance (therefore i?AndNE = 6.6 kpc, Dmzi = 43 kpc), (M/L)AndNE ~ 30(M/L)m3i- This is 
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high enough to argue that, if And NE is at the same distance as M31, it is likely undergoing 
tidal disruption. Yet, And NE may be rather more distant from M31: if instead we assume 
that And NE is ±0.15 mag more or less distant (the error in distance determination), one 
obtains (M/L)A n dNE ~ 6(M/L)m3i- While this latter value still suggests that And NE is 
being tidally disrupted, accurate velocity and distance data will be required to firmly resolve 
the question of whether And NE is bound or unbound. 

4. Conclusions 

SDSS has revealed an overdensity of red giant stars ~ 3° to the NE of the Andromeda 
galaxy; Andromeda NE is within 0^15 in distance modulus (~ 50 kpc) of M31 and is thus 
not an unrelated foreground or background object. Its extrapolated total (/-band magnitude 
is ~ 12.8 (M g = —11.6 or 5 x 10 6 L Q at 760 kpc), and its radius, although not well-defined, 
is ~ 30' (6 kpc). And NE is very diffuse and is probably undergoing tidal disruption, based 
on its distance from M31 and its morphology. And NE possesses a "blue" RGB, unlike 
M31's northern spur, the giant stream, or the outer parts of the present-day M31 disk; yet 
And NE's stellar population is similar to the Gl clump, suggesting perhaps a link with this 
structure. It is possible that it is part of an ancient diffuse stellar stream, or perhaps material 
torn from M31's disk in the distant past. Spectroscopic and deep, high spatial-resolution 
data will be required to further elucidate the nature of this faint, enigmatic stellar feature. 
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Fig. 1. — M31 Stellar Substructure in the SDSS Data: a) The spatial distribution of all stars 
with 21 < % < 22, with an inset optical image of M31 to scale. The arrows indicate major 
halo substructures shown in Ferguson et al. (2002): the northern spur (NS), the giant stream 
(GS) and the Gl clump (G1C). b) Same as a), but with the grayscale adjusted to emphasize 
faint details. The arrow shows the location of Andromeda NE. The orientation of the scan 
is as indicated, c) Same as a), without the inset image, indicating the region used for the 
M31 sample (between the ellipses) and that used for And NE (the solid circle); the dotted 
circle shows the control field used in our analysis. The data have been binned 2' x 2'. X 
and Y are in arcdegrees from the center of M31 (00 h 42 m 44 ! :3 + 41°16'08"5 J2000), along the 
major and minor axes, respectively. The inset image of M31 is from Bill Schoening, Vanessa 
Harvey/REU program/NOAO/AURA/NSF. 
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Fig. 2. — Hess diagrams of M31 and And NE: Top left: Hess diagram of all stars within 
an inner halo ellipse of M31. Middle left: Hess diagram of all stars within the control field 
(see Fig. lc). Bottom left: The difference of the M31 Hess diagram and the scaled control 
field Hess diagram, divided by the square root of the sum of the Hess diagrams. Top right: 
Hess diagram of all stars within a 0?5 radius circle centered on And NE. Middle right: Hess 
diagram of all stars within the control field. Bottom right: The difference of the And NE 
Hess diagram and the scaled control field Hess diagram, divided by the square root of the 
sum of the Hess diagrams. The data were dereddened and transformed to V, I, binned by 
0.2 mag in / and 0.1 mag in (V — I), and are shown with logarithmic grayscales. 
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Fig. 3. — RGB Colors of Stellar Structures from SDSS Data: Left: Hess diagram of And 
NE, minus an appropriately scaled control field Hess diagram, divided by the square root of 
the sum of the two Hess diagrams (to estimate S/N). Left middle: Same as (left), but for a 
field in the Giant Stream (Ibata et al. 2001). Right middle: Same as (left), but for a field in 
the Northern Spur. Right: Same as (left), but for a field in the Gl Clump. The data were 
dereddened and transformed to V, I magnitudes, binned by 0.1 mag in / and V — I, and 
smoothed with a Gaussian filter. The grayscale bar at the top of each panel indicates the 
logarithmic range in estimated S/N shown for that Hess diagram. Fiducial sequences are 
overplotted for Galactic globular clusters with metallicities of (left to right) [Fe/H] = —2.2 
(M15), -1.6 (M2), -0.7 (47 Tuc), and -0.3 (NGC 6553) (Da Costa & Armandroff 1990; 
Sagar et al. 1999), shifted to the adopted M31 distance modulus of 24.4; for NGC 6553, 
(m - M) = 13.7 and E(V - I) = 0.95 were assumed (Sagar et al. 1999). 



